Nuclear spin qubits in a trapped-ion quantum computer 
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Physical systems must fulfill a number of conditions to qualify as useful quantum bits (qubits) 
for quantum information processing, including ease of manipulation, long decoherence times, and 
high fidelity readout operations. Since these conditions are hard to satisfy with a single system, 
it may be necessary to combine different degrees of freedom. Here we discuss a possible system, 
based on electronic and nuclear spin degrees of freedom in trapped ions. The nuclear spin yields 
long decoherence times, while the electronic spin, in a magnetic field gradient, provides efficient 
manipulation, and the optical transitions of the ions assure a selective and efficient initialization 
and readout. 

PACS numbers: 03.67.Lx, 42.50.Dv 



I. INTRODUCTION 

Trapped ions have been recognized as one of the 
most promising candidates for future quantum comput- 
ing (QC) devices for some time [l|. This potential has 
been verified by initial demonstrations implementing sim- 
ple quantum algorithms with trapped ions (see, e.g., H) 
or by entangling eight ultracold ions in a linear trap 'S*!. 
Nevertheless, a number of difficulties have to be over- 
come to scale these devices to larger numbers of qubits. 
The proposals for overcoming these obstacles include seg- 
mented traps where ions are transferred between different 
parts of the trap Q or by joining many spatially sepa- 
rated traps in a quantum network Q . 

Different groups rely on different ions for their experi- 
ments. The NIST group mainly uses ^Be+ @, the groups 
in Innsbruck, Oxford, Atlanta and Wuhan employ ^°Ca+ 
0) H) I9I [lOl . and the group in Maryland works with 
"iCd+ [nf and i"Yb+ [13]. There are also some pro- 
posals and experimental reports using ^''^Yb+ [isl. fl3|. 
43Ca+ [ililll, i38Ba+ [l3l, and In all these 

cases, however, the qubits are encoded in two electronic 
states and the exchange of quantum information between 
ions occurs via quantized vibrational modes. The main 
conditions for these schemes to work include long life- 
times of the computational states and well resolved vi- 
brational modes. Recent experiments have demonstrated 
that quantum information can be conserved almost per- 
fectly even when the ions are moved [Toj . 

In the context of solid-state QC, several proposals have 
been put forward that suggest the use of nuclear spins for 
storing quantum information 20J . Nuclear spins are very 
well isolated from other degrees of freedom and there- 
fore promise long storage times for quantum informa- 
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tion. On the other hand, the weak interaction with their 
environment also makes initialization and manipulation 
of nuclear spins difficult and leads to low detection sen- 
sitivity. It may therefore be necessary to combine nu- 
clear spins with other degrees of freedom to combine 
long storage times with efficient manipulation and de- 
tection [IJ [13, mHi, [H]. Of these proposals, so far 
only one system has been demonstrated experimentally, 
which combines the electronic spin of the diamond NV 
center with a neighboring ^'^C nuclear spin [2^. [27|. 

In this paper, we propose to combine nuclear spins 
with trapped ions. The goal is to design a system that 
combines long decoherence times with efficient manip- 
ulation and detection. The trapped-ion system that we 
consider is closely related to the proposal by Mintert et al 
dH, which rehes on a magnetic field gradient to address 
different ions and to induce couplings between the ions 
that do not rely on the motional degrees of freedom. This 
scheme can be used with different ions, such as ''^Ca^, 
whose nuclear spin is / = 7/2, ^^'''Ba^ and ^^^Ba^ with 
1=3/2 as well as ^^^Yb"'' with / — 1/2. For specific dis- 
cussions, we will refer to the '^'^Ca^ system, but all our 
results are also applicable to the other three candidates. 

The hyperfine splitting of the S'1/2 ground state of 
^^Ca"*" has been precisely measured [2^ js^. Since our 
scheme works best when nuclear and electronic spins can 
be clearly distinguished, we explicitly consider the case of 
a strong magnetic field (i.e., the Paschen-Back regime). 

The paper is structured as follows: In the following sec- 
tion, we first consider a single trapped ion and show how 
the quantum information can be encoded in the nuclear 
spin degree of freedom and transferred to the electronic 
spin. Section III introduces two-qubit operations and in 
section IV, we discuss the operation of the full quantum 
register, including addressing of individual qubits. Sec- 
tion V presents a short summary and conclusions. 
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II. SINGLE *^CA+ ION 

We first consider a single trapped ■*^Ca+ ion. In the 
Paschen-Back regime, we can write the Hamiltonian as 



ESR Transitions 



Ho = nsS, - ilih + AS J, 



(1) 



where we use units of ft = 1. Sz and are the spin 
operators of the electron spin (5 = 1/2) and the nuclear 
spin (/ = 7/2). fig = gsHsB and fli = gn^J^BB are 
the Larmor frequencies of the electron and nuclear spins, 
and A = —806.4 MHz 30] is the hyperfine coupling con- 
stant. We include here only the high- field truncated part 
of the hyperfine coupling. In the absence of external ra- 
diation, the vibrational mode is not coupled to the spin 
degrees of freedom and was therefore not included in the 
Hamiltonian of Eq. (1). 




FIG. 1: Angular momentum states in the ground state of a 
^^Ca"'' ion in a magnetic field B = 1 T (i.e. in the Paschen- 
Back regime) . The four levels drawn in bold are used to define 
two qubits. 

Figure [T] shows the resulting energy level structure in 
a magnetic field of _B = 1 T. The electron Zeeman in- 
teraction causes a splitting of 28 GHz. Each of these 
electronic states is split into 8 nuclear spin sublevels by 
the hyperfine interaction. Due to the nuclear Zeeman in- 
teraction, the splittings within the two multiplets differ 
by 5.7 MHz. 

Transitions between these spin states can be induced 
by the magnetic dipole interaction. Changes of the elec- 
tron spin correspond to frequencies in the microwave 
range, centered around 28 GHz, as shown in the upper 
part of Fig. [21 The lower part shows the transitions cor- 
responding to changes of the nuclear spin state. Their 
frequencies correspond to \A/2\ ± il/. 

While this system could in principle store 4 qubits (i.e., 
three in the nuclear spin and one in the electronic spin), 
we will consider here only two qubits. The first qubit 
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FIG. 2: Magnetic dipole transitions in the *''Ca"'" spin sys- 
tem. The upper trace shows the transitions that correspond 
to a change of the electron spin. They fall into the microwave 
range, close to 28 GHz, and are split by the hyperfine in- 
teraction into 8 transitions, each corresponding to a single 
nuclear spin state. The lower trace shows the transitions of 
the nuclear spin, which fall into the radio frequency range. 
The lower frequency line corresponds to the ms=-l/2 state 
of the electron spin, and the higher frequency transition to 
the ms=+l/2 state. 



corresponds to the electron spin and the second qubit 
will be encoded in the m/ = ±7/2 states of the nuclear 
spin. 

To apply gate operations to these qubits, we use reso- 
nant microwave radiation. 



Hi= H(3 + 2 COs{Vln-,^t)[uJsSx + LOlIx] 



(2) 



Here, Hq is the system Hamiltonian of Eq. ([T]), los and 
LJi are the Rabi frequencies of the electronic and nuclear 
spins and fimw the frequency of the applied field. 

In the following, we will assume that the microwave 
field is weak enough that it interacts predominantly with 
a single transition. This is fulfilled if the Rabi frequen- 
cies are small compared to the splittings between the 
resonance lines, 

ujs ^ 1^1 and ^ . 

The electron spin qubit can then be manipulated by mi- 
crowave radiation, if the microwave field is resonant with 
one of the two transitions at 25.2 and 30.8 GHz, which 
correspond to the nuclear spin being in the to/ = 7/2 
and -7/2 states, respectively. A microwave pulse at one 
of these frequencies thus corresponds to a single qubit ro- 
tation of the S'-spin, conditioned on the /-spin being in 
the corresponding state. Unconditional single-qubit op- 
erations on the S'-spin can be accomplished as a sequence 
of two (or two simultaneous) conditional operations. 

Applying single qubit gates to the nuclear spin qubit 
is less straightforward: In general, a resonant RF pulse 
takes the nuclear spin out of the state space associated 
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with the nuclear spin qubit (to/ = ±7/2). An important 
exception is the operation Pj — e^^'^^^ , which simply 
interchanges the two states |7/2) and | — 7/2) : 

If we use resonant RF irradiation, the exchange be- 
comes again conditioned on the state of the electronic 
qubit. As an example, a RF pulse at Qmw = 406.1 MHz 
exchanges |7/2)^ |l/2)g <^ |-7/2)^ |l/2)g, while the 
\7/2)j \ — l/2)g and |— 7/2)^ | — 1/2)^ states remain invari- 
ant. This corresponds to a CNOT operation, which we 
call CNOTgi: the first index refers to the control qubit 
and the second to the target qubit. Similarly, a reso- 
nant TT pulse applied to the ESR transition at 25.2 GHz 
corresponds to a CNOTig operation. 

We use these operations to implement a SWAP oper- 
ation between electronic and nuclear qubits: 

SWAP = CNOTis CNOTgi CNOTis 
= CNOTgi CNOTis CNOTsi . 

The SWAP operation can now also be used to imple- 
ment arbitrary single-qubit operations Uj on the nuclear 
spin qubit: 

Ui = SWAP Us SWAP , 

where Us is the single-qubit operation acting on the elec- 
tron spin qubit. The system allows thus to implement 
arbitrary one- and two-qubit operations. 

III. A PAIR OF QUBITS 
A. System and Hamiltonian 

To implement a quantum register on the basis of the 
"'^Ca^ system, we need a string of trapped ions. We will 
use the | ±7/2) states of the nuclear spins as the qubits, 
while the electron spins are used to couple the nuclear 
spin qubits to each other. 

When a magnetic field gradient B{z) ^ Bq + b6z is 
applied with b — dB/dz, the Hamiltonian of the system 
is 

H2 = Y.n'sSi-Y,n',ii + Aj2siii 

i i i 

-J2j,,SlSi, (4) 

where the first and the second terms represent the Zee- 
man energies of the electronic and the nuclear spins, re- 
spectively, and the third term is the hyperfine interac- 
tion. The Larmor frequencies fi^ and depend on the 
magnetic field gradient b and the position Zi of the ions 
as 



and correspondingly for the nuclear spins. 

The fourth term in Eq. ^ represents a coupling be- 
tween the different ions, which takes the form of an Ising 
interaction. This effective interaction is generated by the 
magnetic field gradient. The coupling constants are 
proportional to the square of the ratio of the magnetic 
field gradient to the center of mass frequency [H, [3l|, Is^l ■ 
It determines the gate operation time T = tt/ J12 for two- 
qubit operations. In TableU we list some possible combi- 
nations of trap parameters together with the correspond- 
ing gate operation times. For simplicity, we have omitted 
the corresponding interaction between the nuclear spins, 
which is very small compared to all other terms. 

We first consider the case of two ions in the trap. For 
the parameters corresponding to the first line in Table 
m the spectrum of the ion pair remains essentially the 
same as that of the single ion, as shown on the bottom 
of Fig. [21 The separation of the resonance lines due to 
the magnetic field gradient becomes apparent on an ex- 
panded frequency scale, as shown in the central trace 
of Fig. [31 The splitting due to the gradient-induced J- 
coupling between the ions is then another 3 orders of 
magnitude smaller, as shown in the top trace. 

TABLE I: CNOTsiSa gating time T for a pair of trapped 
''^Ca^ in a 1 T magnetic field for different center-of-mass 
trap frequencies vi, magnetic field gradients b, and distances 



vi/2n (MHz) 


b (T/m) 


AZmin i^J.m) 


J12 (kHz) 


T (ms) 


1.0 


450 


5.5 


3.25 


0.24 


1.0 


230 


5.5 


0.85 


0.92 


1.0 


50 


5.5 


0.40 


1.75 


0.8 


340 


6.3 


2.90 


0.27 


0.8 


160 


6.3 


0.65 


1.21 


0.8 


35 


6.3 


0.30 


2.62 



B. Two-qubit operation 

A CNOT operation between the two qubits can be 
achieved via the following three steps: 

(1) The information encoded in nuclear spins is trans- 
ferred to the electronic spins by the gates SWAPiiSi and 
SWAPi.s^ . 

(2) CN0Tgjg2 is apphed to the electronic spins. 

(3) Step (1) is repeated to transfer the information 
back to the nuclear spins. 

Looking in more detail at the second step, we note that 
the CNOTg^g^ operation does not involve the nuclear 
spin. We may thus disregard the nuclear spin degrees 
of freedom during this step. The problem reduces then 
to the conventional system of two spins 1/2, which are 
described by an effective Hamiltonian 



fls = giJ-BiBo + bzi) 



(5) 
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Here, the Larmor frequencies ^1.2 of the two electron 
spins include a contribution from the hyperfine coupling, 

and analogously for the second ion. Here, m/^ = ±7/2 
denotes the nuclear spin state. This Hamiltonian is iden- 
tical to that of a weakly coupled pair of spins 1/2 in 



J 



liquid-state NMR. The CNOT operation can therefore 
be carried out in complete analogy, either by a selective 
pulse on one of the resonance lines correspon ding to the 
transition between the states (1,0) and (0,1) [3311, or by 
two tt/2 pulses applied to the first qubit, separated by a 
free precession period of duration n/ J12 under the cou- 
pling Hamiltonian [s^, HE] • 




Frequency (GHz) 



FIG. 3: ESR spectrum of two ions in a trap. The lowest trace 
shows the full spectrum, the middle trace shows the lowest 
frequency line on an expanded scale. The splitting is due to 
the magnetic field gradient. The uppermost trace shows the 
splitting due to the J-coupling between the ions. 



C. Selective SWAP operations 



D. Initialization and readout 



The SWAP operations required in steps (1) and (3) 
were discussed in section |TT1 However, since the quan- 
tum register now contains multiple ions, we need more 
specific control over the radiofrequency and microwave 
fields. Since the overall duration of any algorithm is 
determined mostly by the relatively weak J-couplings, 
it is not necessary to make the Rabi frequencies very 
high. A useful Rabi frequency for the microwave pulses 
would then be of the order of 1 MHz, which can easily 
be achieved with very small microwave power. As an ex- 
ample, the devices described in [3^ can generate 1 MHz 
Rabi frequencies with 6 fiW of microwave power. Being 
much smaller than the hyperfine coupling constant, the 
microwave will excite only the transition corresponding 
to the desired nuclear spin state and only the targeted 
ion, while the neighboring ions will not be affected. Such 
a pulse will, however, excite both lines of the doublet due 
to the J-coupling. 



Apart from the gate operations discussed so far, ini- 
tialization and readout are essential steps for the imple- 
mentation of a quantum information processor. In this 
respect, our proposal is almost identical to that of con- 
ventional ion trap quantum computers (37j : Initialization 
is achieved by optical pumping, while readout is achieved 
by electronic shelving amplification [38]. Since the ad- 
dressing of the ions as well as the interactions between 
them are magnetic, there is no need for cooling the ions 
to their ground state [l^, but only into the Lamb-Dicke 
limit. 



The main difference arises from the fact that the qubits 
are encoded in the nuclear spin degrees of freedom: After 
the initialization, it is thus necessary to SWAP the ini- 
tialized state into the nuclear spin. The same process is 
required in reverse for readout: The resulting state must 
be SWAPed from the nuclear spin degree of freedom to 
the electronic degree of freedom, where readout occurs. 
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IV. QUANTUM REGISTER 
A. Addressing qubits 

For an actual QIP implementation, more than two ions 
are required. In terms of trap dynamics, it should be pos- 
sible to store more than ten ions in such a trap. Address- 
ing of the ions is straightforward in the scheme described 
here, as long as the Larmor frequencies of the electronic 
spins do not overlap. This requires 

g^j,BbAzN<A, (6) 

where N is the number of ions in the string. 

Addressing the individual qubits can then be achieved 
by keeping the microwave field strength well below the 
difference of the Larmor frequencies of two neighboring 
ions, 

uji < g fisb Az . (7) 

This condition can be easily met without increasing the 
total duration of any quantum algorithm. 

So far, we have only considered the electronic spins. 
But as we have discussed in the previous section, actual 
gate operations require resonant excitation of the elec- 
tronic as well as the nuclear spin. The magnetic field 
gradient also separates the nuclear Larmor frequencies 

by 

Ani = -mbAz, (8) 

with Qj in MHz, b in T/m, and Az in m. Compared to 
the electron spin resonance frequencies, this splitting is 
reduced by approximately 4 orders of magnitude, which 
amounts to « 700 Hz for the parameters used above. 
Addressing the nuclear spins individually would thus in- 
crease the duration of a quantum algorithm significantly. 
However, this can be easily avoided if we notice that only 
TT pulses are applied to the nuclear spins and that all rel- 
evant operations discussed in sections II and III involve 
two SWAP operations between the electron and nuclear 
spins of the ion and thus an even number of it pulses 
applied to the nuclear spins. We can therefore apply the 
RF pulses nonselectively, i.e. to all nuclear spins in the 
quantum register, using a strong rf field and short pulse 
durations. The individual vr pulses will then affect all nu- 
clear spins in the trap; however, for those ions that are 
not addressed by the frequency-selective pulses applied 
to the electronic spins, the operation of the tt pulse is 
reversed by the next tt pulse and results in a NOP for 
the inactive qubits. 

B. Refocusing operations 

In a string of ions, the induced coupling varies with the 
position of the ions in the trap [Tsl . [28j . Couplings exist 
between all pairs of ions in the trap, not only between 



nearest neighbors. Quantum algorithms therefore have 
to take the full network of couplings into account. This 
operational form of the couphngs, SIS^, provides a close 
analogy to liquid-state NMR, where undesired couplings 
are usually eliminated by refocusing techniques (see, e.g., 
[3l,|3l[3i,|40,lill,|4l). 

In the present case, such refocusing operations are 
more straightforward than in the case of liquid state 
NMR. If we want to implement a two-qubit operation be- 
tween ions i and j, we want to have one active coupling 
between these two ions, while all other couplings should 
be refocused. Since all terms in the relevant Hamilto- 
nian in Eq. (j4|) commute with each other, it is sufficient 
to consider the terms 

H^^n'sSl+^n,Si-J,,SlSi-A,SlS^, (9) 

where the index k represents any of the passive spins. 
The evolution of the two active spins i and j under the 
coupling Jij is now achieved by letting the system evolve 
under this Hamiltonian for a duration r = 7r/(4Jij), ap- 
plying TT pulses to spins i and j, and letting them evolve 
for another period r, and applying another pair of tt 
pulses. The resulting overall evolution is then 

[/2 ^ e-»-ff4(^/4J.,)g-»7r(S;,+S^)g-»_f/4(^/4J.j)g-i7r(Sj,+SJ) 

(10) 

Using 

~^S^z ~ ^S^z ~ Jij^z^z + JikSlS^ , (11) 

we can simplify this to the desired evolution 

U2 = e^'^-^^^'- . (12) 

C. Decoherence 

Like in other implementations of quantum information 
processing, the operations must be completed in a time 
short compared to the decoherence time. We consider 
here two contributions to the decoherence: (i) dephasing 
due to magnetic field instability and (ii) heating of the 
motional degrees of freedom. 

The evolution of the electronic spins in the magnetic 
field changes the phase of the qubits by 

Aipt = QiT . 

Instabilities of the magnetic field therefore affect the 
phase of the overall quantum state and must be reduced 
to 

AB « . 
gi^BT 

Since the gate operation time r is of the order of a few ms, 
this requires a stability of the order of 10~* T. Clearly, 
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this is a challenging condition; however, the correspond- 
ing magnet technology is well established in the field of 
high-resolution NMR, where field stability of better than 
10"^*^ is required. The requirements can also be reduced, 
by working at lower magnetic field, by refocusing tech- 
niques, or by working in decoherence-free subspaces [isl ]. 

Additional contributions to the decoherence of trapped 
ions are due to mobile charges in the trap and to mo- 
tional heating. In the proposed trap, where the infor- 
mation is encoded in the spin degrees of freedom, such 
contributions should be considerably smaller than in con- 
ventional ion traps. Moreover, the couplings between dif- 
ferent electron spins are due to the Coulomb interaction 
of the trapped ions under the magnetic field gradient, 
which only virtually excites the motional degrees of free- 
dom. This also implies that motional heating during gate 
operation is suppressed, as discussed by Wunderlich and 
coworkers [l^, [2^ . 



D. Experimental feasibility 

A main point of our proposal is that we use only spin 
degrees of freedom. This leads to a considerable sim- 
plification of the level diagram compared to the case 
where electronic degrees of freedom contribute to the 
relevant Hilbert space spanned by 6*1/2 8'iid I?5/2 p^ . 
The possibility of addressing trapped ions using a mag- 
netic field gradient and resonant microwave radiation was 
recently demonstrated using a magnetic field gradient 
weaker than 1 T/m [ij, |T1]. Although such a mag- 
netic field gradient is not big enough to generate efficient 
CNOT gates, this should become possible by an opti- 
mization of the experimental parameters. Stronger mag- 
netic field gradients (up to 8000 T/m) can be generated, 
as discussed by Mintert and Wunderlich |28i] . 

A larger field gradient would provide stronger cou- 
plings and therefore faster two-qubit gate operations. 
However, the gradient also couples the electron spin to 
the vibrational motion and thereby, if it is too big, re- 
duces the fidelity of the gate operations. An upper limit 
to the useable gradient strength is given by the condition 
e <^ 1 [2^ . where the parameter 

in units of ^ 1 is a generalized Lamb-Dicke parameter, 
with m the mass of the ion, N the number of the confined 
ions, and vi the fundamental vibrational frequency. For 
our system, where vi — 2tt MHz, = 2, m 43 x 1.67 x 
10~^^ kg the upper limit for the magnetic field gradient 
becomes h < 455 T/m. Limitations on the magnetic field 
gradient for stable trapping of ions are also discussed by 
Cronin et al. Q. 



V. DISCUSSION AND CONCLUSION 



In this paper, we have proposed a scheme for encoding 
quantum information in nuclear spin states of trapped 
ions. '^■^Ca+ has also been used in Ref. (IQ], where the au- 
thors encoded the information in field-insensitive hyper- 
fine states. Since we work in the Paschen-Back regime, 
there are no field-insensitive states. It is therefore nec- 
essary to use stable magnetic fields and, possibly, shield 
against undesired external fields. However, the require- 
ments on the stability of the magnetic field remain signif- 
icantly less strict than in the case of liquid-state NMR, 
where the field has to be kept stable on a scale of w 10"^". 
We conclude that this should not be a significant obsta- 
cle. 

The proposed scheme is also related to an earlier pro- 
posal for QC with doped fuUerenes [13, [13, ll^- Both use 
nuclear and electronic spin degrees of freedom and work 
in the Paschen-Back regime. In terms of requirements for 
field stability both systems should be comparable and 
the gate operations rely in both cases on resonant mi- 
crowave pulses. The main advantage of the trapped-ion 
system over the fuUerene scheme would be the relatively 
easy readout, which remains a significant challenge for 
the fullerene system. While the solid-state systems may 
be easier to scale up, our trapped-ion system may also be 
made scalable by approaches based on segmented traps 
[1| or multiple connected traps A third possibihty 
might be the use of microtraps for single ions, which is 
also compatible with field-gradient induced couplings be- 
tween the ions fs^ . 

Throughout this paper, we have relied on field-gradient 
induced couplings between the ions to implement multi- 
qubit gate operations. This is not the only possibility, 
however, because the use of nuclear spin degrees of free- 
dom for shelving the quantum information is also com- 
patible with the more conventional Cirac-ZoUer method 
in or the S0rensen-M0lmer scheme [i^ which does not 
require extensive cooling of the ions. 

In summary, we have investigated the possibility to 
carry out QC using nuclear-spin qubits in '*'^Ca+ ions. 
Using the nuclear spin should result in long decoherence 
times and working in the Paschen-Back regime results 
in a straightforward distinction of nuclear and electronic 
spin degrees of freedom. All the required gate operations 
can be achieved with microwave pulses, in which indi- 
vidual addressing of the qubits is achieved by frequency- 
selection. If a few-qubit quantum computer on the basis 
of this scheme can be implemented, it could serve as a 
test-bed not only for larger trapped-ion quantum com- 
puters, but also for similar schemes based on condensed- 
matter systems, such as endohedral fuUerenes p^[23l.[25j 
or NV centers in diamond [1^, [l^l . 
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